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Abstract 

Zone-center phonon frequencies of polar lattices are calculated for uniaxial 
crystals from the symmetry arguments. Long-range Coulomb forces and crys- 
tal anisotropy are explicitly taken into account. Free-carrier contributions into 
a dielectric constant are included. The angular dispersion of optical-phonon 
modes is compared with data for hexagonal 6H-SiC politype. 
PACS: 63.20.-e, 78.30.-j 

I. INTRODUCTION 

Electrostatic dipole-dipole interactions play an important role in the theory of lattice 
vibrations. It is common knowledge |l[ that the degeneracy of phonon modes at the Brillouin 
zone-center (e.g., in the cubic 3C-SiC crystal) is removed if the atomic displacements are 
accompanied by the Coulomb field. Then the frequency of the longitudinal optical mode 
becomes larger than the frequencies of transverse modes. For noncubic crystals (e.g., for the 
hexagonal or rhombohedral SiC polytypes), the long-range Coulomb field gives rise also to 
an angular dependence of the zone-center modes: at k = the optical-phonon frequencies 
depend on the propagation direction. 

Such a phenomenon is rather unusual from both the physical and mathematical point of 
view: the eigenvalues of dynamical matrix calculated for k = depend on the k-direction. 
This is caused by the nonanalytic k-dependence of the dynamical matrix which results from 
the long-range dipole-dipole interaction. In polar cubic crystals, the Coulomb field splits 
the three-fold degeneracy of optical modes at the Brillouin zone-center, but the frequency 
dependence on the propagation direction as well appears in uniaxial crystals due to the 
long-range electrostatic field. 

The electro dynamic part of the problem was formulated by Loudon M. The Coulomb 
contributions in the dynamical matrix are usually calculated be means of an Evald sum- 
mation The angular dispersion of optical modes is clearly demonstrated by the recent 
numerical calculations of the zone-center phonons || and for the whole Brillouin zone 0] 
in the case of semiconductors A 3 B 5 with the wurtzite structure. The Coulomb field is also 
taken into account in the theory of phonon-plasmon coupled modes (polaritons) ||, when 
the effect of free carriers is studied. 
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The main aim of this paper is i) to calculate the angular dispersion for the zone-center 
phonons in uniaxial crystals using the symmetry arguments and ii) to consider the effect of 
free carriers on these modes. For the sake of definiteness, we are interested in the phonon 
modes of uniaxial SiC polytypes which are very popular now in technical applications. 



II. OPTICAL MODES AT THE ZONE-CENTER OF CUBIC CRYSTALS 

Among the hexagonal and rhombohedral SiC polytypes, there is the cubic 3C-SiC one 
with two atoms in the unit cell. First we consider the optical modes in this simplest case. For 
the nearest vicinity of Brillouin zone-center k <C ir/d, where d is the lattice parameter, the 
acoustic and optical modes can be divided using the series expansion in k of the dynamical 
matrix. As the result, we obtain in the zero approximation in k for the optical displacements 
Ui {i = x, y, z) the system of three equations 

(0 - M*w 2 )u = f, (1) 

where M* is the reduced mass of two atoms (Si and C) in the unit cell, <fi is the diagonal 
element of the force-constant matrix (the only one diagonal element of the 3 x 3-matrix exist 
in a cubic crystal). The value of <j) can be calculated in the nearest-neighbor approximation. 
As for the long-range Coulomb interaction, it can not be considered in such a way. The 
Coulomb effect is described by a force f = ZeE acting on an effective charge Z, where 
the electric field E is found from Maxwell's equations. Eliminating the magnetic field from 
Maxwell's equations, one can express the electric field E in terms of polarization P: 

-4vr[k(kP)-^P/c 2 ] 

k 2 -LU 2 /c 2 ' 1 ' 

We are interested in uj of the order of optical mode frequencies, i.e., uj/c~ 10 3 cm -1 . If 
the phonon is exited by light, its wave vector has the value of the photon wave vector, i.e., 
of the order of 10 5 cm -1 . Then the condition k ^> u/c is fulfilled, and the terms with c 2 
have to be omitted in Eq. (Q) which becomes: 

E = -47rk(kP)/fc 2 . (3) 

In the long- wave limit (k <C n/d), the polarization is related to the phonon displacement 
and electric field by the macroscopic equation: 

P = NZeu + xE (4) 

where \ is the atomic permittivity and N is the number of unit cells in 1 cm 3 . Sometimes, 
the local field is used in equations like (f|) instead of the macroscopic field E. But for cubic 
crystals (when the simple Lorentz relationship exists only), the local field can be eliminated 
renorming the force constant (p. 

Using Eqs. (|3[) and (|j), the electric field E can be expressed in terms of u. Then Eq. 
(0) gives the frequencies of transverse and longitudinal optical modes in the cubic crystal 

u 2 TO = <p/M* and u 2 LO = <p/M* + p (5) 
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where 



p = Z 2 e 2 N/e°° and e°° = 1 + 4ir X - 



(6) 



Although the relation 
quencies of optical modes 
a cubic crystal. 



between E and P involves the k-direction explicitly, the ire- 
are independent of the propagation direction as it must be for 



III. OPTICAL MODES AT THE ZONE-CENTER OF UNIAXIAL CRYSTALS 



The crystal anisotropy of the noncubic SiC polytypes is known to be small because the 
nearest neighbors of any given atom conserve the cubic symmetry. Let us introduce the 
strain tensor which describes a small difference between the dynamic matrices for the 
noncubic polytype and the cubic one. Then the phonon spectrum of the noncubic polytype 
can be obtained in the following way. In the first step, we transform the Brillouin zone of 
the cubic polytype ("the large zone") using the strain e^. Hence, we find the frequencies of 
the so-called strong modes. For the zone-center, they can be obtained by the expansion of 
the dynamic matrix in the strain e^. 

In the second step, we take into account that noncubic polytypes have more than two 
atoms in the unit cell and the additional optic modes appear. Phonon branches of the large 
zone are folded || into the Brillouin zone of the noncubic polytype giving additional weak 
modes. The weak-mode intensity in optics and Raman scattering was calculated in 
Therefore in the present paper, we consider only strong modes. 

The dynamic matrix can contain only components invariant under the symmetry 
transformations of the crystal. There are two invariants of first order, e zz and e xx + e yy , the 
z-axis being parallel to the c-axis. We can fix the crystal volume, i.e., impose the condition 
en = 0. Then we have only one invariant, for instance, e zz which can be written only on the 
diagonal of the force-constant matrix of Eq. (|I|). The coefficients of the xx and yy elements 
are equal because of the rotation invariance around the c-axis. At last, we can omit the 
common frequency shift. Therefore instead of Eq. ([I]), we obtain 



ft + pn 2 x - u 2 
pn x n y 
pn x n z 



pn x n y 
ft + pn 2 - uo 2 



pn y 



n, 



pn x n z 
pn y n z 
a + pn 2 — oo 2 



( u x 

Uy 

V u 2 



0. 



(?) 



where n = k/fc and 



a 



<f)/M*, ft = a + be z 



(8) 



We take the vector k in the ?/z-plane and denote 9 the angle between k and the c-axis: 
n x = 0,n z = cos 8 and n y = sin 8. Then we see from Eq. (0) that there are one transverse 
mode (TOi) vibrating in x-direction and two modes in the yz-plane with the frequencies 
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y.z 



(0) 



Up + a + ft) ±)-{[p + {a - ft) cos2#] 2 + (a - ft) 2 sin 2 2d} 



TOI 

2 : 



ft, 

1/2 



(9) 
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Emphasize that Eq. @ gives the phonon frequencies in the zone-center, but they depend 
on the propagation direction 9. This dependence has its origin in the simultaneous effect of 
the Coulomb field (described by the constant p) and crystal anisotropy {(3 ^ a). In absence 
of the Coulomb field (p = 0), we have u 2 = a, u 2 = (3, and there is no angular dispersion. 
For the isotropic case (a = 0) , Eq. (|9|) gives the modes for the cubic crystal. 

If the Coulomb effect is small in comparison with the crystal anisotropy (p \a — /3\), 
we can omit the off-diagonal terms in the matrix ([?]). So there are one mode vibrating close 
to the c-direction with the frequency uj 2 = a + pcos 2 9 (with an accuracy to p 2 /(a — /3) 2 ), 
and the other mode near the ^-direction with the frequency u 2 = (3 + psin 2 9. 

In the opposite limiting case of the small crystal anisotropy, it is useful to pass to the 
coordinate system with the z'-axis along the k-vector, making in Eq. fl7|) the unitary trans- 
formation 

/I \ 
Uij = cos9 sin^ (10) 
\ — sin 9 cos 9 J 

Then we arrive to the diagonalization problem of the matrix 

/ (3 \ 

(3 cos 2 9 + a sin 2 9 ((3 — a) sin9 cos# 
\ ((3 - a) sin 9 cos 9 (3 sin 2 9 + a cos 2 9 + p J 

We see that besides the TOi mode, in the case \a — /3| p, there are another nearly 
transverse TO2 mode and the nearly longitudinal LO mode with the frequencies 

u4 O2 (0) = Pcos 2 9 + asm 2 9, (12) 
uj 2 lo (9) = p + (3 sin 2 9 + a cos 2 9, 

which can be obtained also by the expansion of Eq. (||) with an accuracy to (at — (3) 2 / p 2 . 
The dispersion curves (0), ( |i~2"D are shown schematically in Fig. 1. The angular dispersions 
of the form (^) were obtained by Loudon 0. 

One can see from Eq. ([5|) that a conservation law exists. Namely, the sum of the squared 
frequencies of the yz modes is independent of the propagation direction, e.g., 

^ = 0) + uj 2 (9 = 0)=u 2 y (9 = tt/2) + u 2 z (9 = vr/2). (13) 

As an example let us consider 6H-SiC polytipe. The angular dispersion of its optical modes 
is known from the experiment U, @. For 9 = (propagation parallel to the c-axis) the TOi 
and y mode are degenerate, and their frequencies are equal to y/]3. The experimental value 
is 797 cm -1 (with uncertainty about 1 cm" 1 ). The corresponding value for the longitudinal 
mode is ojlo{9 = 0) = y/p + a. For 9 = n/2 (propagation perpendicular to the c-axis) 
Ut02(9 = 7r/2) = y/a (experimental value is 788 cm -1 ) and ljlo(0 = vr/2) = y/p + [3 (970 
cm -1 ). It immediately follows that: p = 552. 9 2 cm -2 , a = 788 2 cm -2 and (3=797 2 cm -2 . 
Thus one calculates ujlo{6 = 0) = y/p + a =962.6 cm -1 , which should be compared with 
the experimental value 964 cm -1 . The small difference between these two values can be 
attributed to the anisotropy in the atomic permittivity which is considered in the following 
section. 

4 



(11) 



IV. EFFECTS OF THE PERMITTIVITY ANISOTROPY AND FREE CARRIERS 



We assumed in the previous section that the uniaxial anisotropy affects only the short- 
range contribution to the force-constant matrix. But in uniaxial crystals the atomic per- 
mittivity x is a tensor with two independent components x\\ an d X± corresponding to the 
crystal axes. This effect is small because each atom has nearly cubic surroundings, but it 
should be included for a careful comparison with experiments. In a similar way, free carriers 
make a contribution to the angular dispersion of the longitudinal optical mode. 

To take into account both the anisotropy of atomic permittivity and the conductivity of 
free carriers a, we write instead of Eq. (PI) the following 



(14) 



P± = NZeu ± + ( X ± + E ± . 



Using Eqs. (0) and QI4"), we obtain the equation of motion in the form (0) and the 
phonon frequencies @ , but the conservation law (|13|) does not work now because p becomes 
the function of 6: 



p{6) = Z 2 e 2 N 



.Cll 



OJ 



ef + Am^ ) cos^ 6 + [ ef + 4vrz 



.0"_L 



OJ 



sin 2 6 



(15) 



where = 1 + 4ttx||, £j° = 1 + 4:irx±- Notice, that the vibration modes obtain same 
damping due to conductivity. In addition, the optical phonon has a natural width F given 
by its probability to decay into phonons of lower energy, and the term iF should be added 
to u in Eq. (0). 

Then we can use transformation ([H]) and obtain matrix (|ll|) with the function p(6) 
instead of constant p. We see that in the limiting case of the weak anisotropy, |a — /3\ 
p(0), the Coulomb field (and therefore the carriers) affects only the longitudinal mode. Its 
frequency is determined by the equation 



R(u) = p(9) + (3 sin 2 9 + a cos 6 — iujT 



UJ 



0. 



(16) 



where p{6) given by Eq. (|T5|) depends on to explicitly and by means of the conductivity a. 

Equation ([RJ) gives the frequency of the LO-phonon-plasmon coupled mode in uniaxial 
semiconductors. Notice that in the isotropic case, Eq. ( |T6"D coincides with the condition 
e(uj) = 0, where the dielectric function e(u) is given by the well-known expression 



1 + 



UJ 



to 



UJ. 



TO 



UJr, 



UJ, 



TO 



— UJ' 



iujT uj{uj + i^) 



and the plasmon frequency 



oj 2 = Anne 2 /e°°m. 



uj 2 to + Z 2 e 2 N/e°°, and the 



In this case, Eqs. (||), (||), and (g) give uj^ = a = j3, uj lo 
Drude formula for conductivity reads a = ne 2 /m(—iuj + 7). 

The function R(uj) from Eq. (^) is measured in Raman experiments. Namely, the 
Raman intensity as a function of frequency transfer oj is 
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(17) 



for the LO mode excitation with the propagation direction 6. If the incident or scattered 
light has a finite aperture, Eq. flTTD should be integrated over allowed 6. 

Eq. ([17]) can be used in experimental studying the effect of carriers on the Raman 
scattering in uniaxial semiconductors. Then the conductivity tensor in Eq. ( |Io] ) is given 
by the Drude-like formula with the diagonal components m^± and 7||,_l, for instance, cry = 
ne 2 /m\\(— iu + 7y). 

Let us emphasize the main result of the paper: the effects of crystal anisotropy (a ^ 0) 
and Coulomb field p(6) on the phonon dispersion are explicitly separated as one can see in 
Eqs. (D, (£§. 
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Figure caption Angular dispersion of optical-phonon modes in uniaxial crystals at the 
zone-center. The angle 9 is the angle between the c-axis and the wave vector k — > 0. The TOi 
mode is polarized perpendicular to the c-k plane. The LO and T0 2 modes have a nearly 
longitudinal and transverse character, respectively, if the Coulomb force effects dominate 
over the crystal anisotropy. 
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